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Mass splittings and matrix elements of mesons and baryons containing a 
single heavy quark 

UKQCD Collaboration - presented by David Richards a 

a Department of Physics and Astronomy, The University of Edinburgh, Edinburgh EH9 3JZ, Scotland 

We present a study in the quenched approximation of the B parameter Bb and the decay constant fs using 
heavy-quark propagators implemented in the static approximation, and light-quark propagators computed using 
an 0(a)-improved fermion action. We find a value of Bg 3 " 1 close unity, and discuss the systematic errors entering 
into the calculation. /|| at is extracted using a variational fitting technique in order to obtain a reliable estimate 
of the ground state. 

In the second part of the talk, we describe an exploratory study of baryons containing a single heavy quark, 
computed using the 0(a)-improved fermion action. We obtain masses generally in good agreement with exper- 
iment in both the charm and beauty sectors. We also report preliminary results for the form factor G\ in the 
semi-leptonic Ab — > A c transition. 



1. INTRODUCTION 

The investigation of hadrons containing a single 
heavy quark affords the opportunity to determine 
the least well known elements of the CKM ma- 
trix. However, their determination is contingent 
upon a precise understanding of the effect of the 
strong interaction on weak decays. This has pro- 
vided the spur to an intense study of heavy-quark 
physics in lattice gauge theory. As well as their 
phenomenological importance, lattice studies of 
heavy-quark systems may also answer important 
theoretical questions such as the nature and ap- 
plicability of the Heavy Quark Effective Theory 
(HQET). In this talk, we present two related lat- 
tice studies of heavy-quark systems. The first is 
an investigation of the leptonic decays and mixing 
parameter of the B system, using heavy quarks 
computed in the static approximation. The sec- 
ond is an investigation of the spectrum of baryons 
containing a single heavy quark, and the first re- 
sults of a study of the semi-leptonic form factors 
of the Aj — ► A c transition. Here we use propa- 
gating heavy quarks implemented using an 0(a)- 
improved fermion action. We conclude this intro- 
duction with an outline of details of the lattice 
simulation common to both studies. 

The calculation is performed on 60 config- 
urations of a 24 3 x 48 lattice at (3 = 6.2. 



The light-quark propagators are computed using 
the Sheikholeslami-Wohlert (SW) tree-level 0(a)- 
improved fermion action [1], with "rotated" prop- 
agators so that the leading discretisation errors 
are 0(aa s ) for light-quark systems [2]. Propaga- 
tors are computed at three values of the light- 
quark mass, corresponding to k; = 0.14144, 
0.14226 and 0.14262. These lie in the region 
of the strange quark, whose mass as determined 
from the ratio m 2 K /m 2 p corresponds to k s = 
0.1419(1); massless quarks correspond to K cr ; t = 
0.14315(2) [3]. 

For dimensionful quantities, we take the inverse 
lattice spacing to be a -1 = 2.9(2) GeV, where 
the quoted error is sufficient to encompass the 
values determined from f w , Ro, m p and the string 
tension. 

2. B B AND /| tat 

The static approximation arises from keeping 
only the leading term in the expansion of the 
heavy-quark propagator [4]. In this calculation, 
the static quark propagators are computed from 
"smeared", or spatially-extended, sources to both 
local and smeared sinks; smearing has proved es- 
sential to obtain a signal in this case [5]. The 
light-quark propagators are computed using local 
sources and sinks. We use both gauge-invariant 



2 



Jacobi smearing [6], with an r.m.s smearing ra- 
dius of ro ~ 6.4, and several Coulomb-gauge-fixed 
smearing functions. 

2.1. Lattice operators and renormalisation 

In the continuum full theory, the pseudoscalar 
decay constant of the B meson is defined through 

{Q\A,(Q)\B(p)) = ip,f B (1) 

while the B-parameter is defined through 

{B°\O L (0)\B°), 
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where Ol is the AB = 2 four-fermi operator 

L = (Ml -75)?) (Ml -75)?), ( 3 ) 
and n is the renormalisation scale. The 
renormalisation-group-invariant (RGI) B param- 
eter is then defined by 

B B =a s {n)- 2 lt } °B B {n). (4) 

The matching of these operators in the full con- 
tinuum theory to the relevant operators in the 
static lattice theory is performed as a two stage 
process: the matching of operators in the full con- 
tinuum theory to operators in the static contin- 
uum theory, followed by the matching of opera- 
tors in the static continuum theory to those in 
the static lattice theory. Here we will illustrate 
the procedure in the case of the B parameter. 

The one-loop matching factor between the con- 
tinuum full theory at a scale mj, and the con- 
tinuum effective theory at a scale /j, < mj 
introduces mixing with the operator Os = 
(6(1 — 75)?) (b(l — 75)?), and assumes the form 
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with C s = -8 and C L = -14 [7]. 

The matching of the static continuum theory 
to the static lattice theory introduces mixing with 
two further operators through the breaking of chi- 
ral symmetry 
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where Dl,Dr and Djv are listed in ref. [7] and 
[8]; note that Dl is corrected to use the reduced 
value of the quark self-energy [9]. Expanding 
Equation (5) and (6) to 0(a s ), and using the 
boosted value of aj, att (a _1 ) [10], we obtain 



Of\m b ) = J2 Z ^ 
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i=L,N,R,S 



with Z L = 0.55, Z N = -0.15, Z R = -0.04 and 
Z s = -0.18. 

We note that the 0(a s ) contributions to the 
matching coefficient are, in general, large, ques- 
tioning the applicability of the one-loop perturba- 
tion theory calculation. Indeed, if we first com- 
pute O e L a (a~ 1 ) according to Equation (6), and 
insert the result into Equation (5), we obtain a 
result for B B some 20 % larger than that quoted 
below, thus providing an estimate of the possi- 
ble uncertainty due to higher order perturbative 
corrections. 

The calculation of Z^ at , the matching factor 
for f B tat , follows the procedure described above, 
with the simplification that there is no operator 
mixing. Here we use Z^ at = 0.78, very close to 
the values employed in our previous study [11] 
and by APE [12]. 

2.2. Calculation of /| tat 



To extract fg 
lators: 



we study the two-point corre- 
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where E is the (unphysical) binding energy, and S 
and L refer to operators computed using smeared 
and local heavy-quark fields respectively. Z L is 
computed from the ratio of C LS and C ss , with 
E and Z s extracted from a fit to C ss . The pseu- 
doscalar decay constant f B tat is then related to 
Z L via 
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In order to ensure that we have isolated the 
ground state reliably, for the case of gauge-fixed 
smearing we compare the values for Zl obtained 
from a singe-state fit, to those obtained to a two- 
state fit to a matrix of correlators [13]. We find 
consistency both between one-state and two-state 
fits, and between fits with differing smearing func- 
tions. We quote as our best estimate the value 
obtained from a two-state fit to the exponentially- 
smeared gauge-fixed correlators, where the statis- 
tical errors are sufficient to encompass the values 
from the different fitting procedures and smearing 
functions. We obtain Z L = 0.112 1 g, yielding 



stat = 26 g 
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where the second error on f B tat arises from the un- 
certainty on the scale. The salient feature is the 
agreement between Zl obtained from simulations 
in which the light quarks are computed using the 
Wilson and using the SW action [12]; the discrep- 
ancy between the central values of f B tat obtained 
with the two actions arises principally from the 
very different perturbative values of Z^ at - 

2.3. B parameter Bb 

We compute the B parameter by evaluating the 
three-point correlators 
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where i = L, R, S, N . The amplitude Zs and 
the exponential factors are eliminated by study- 
ing the ratios 
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We then have 

BT t (m b ) = J2 z ^ z A at r 1 R i 



(14) 
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Figure 1. The B parameter at each k; is shown 
as the squares for gauge-invariant Jacobi smear- 
ing. The solid (dashed) lines show the correlated 
(uncorrelated) chiral extrapolations respectively, 
and the crosses the extrapolated values. 



where the Z{ are given following Equation (7). 

The RGI B parameter at each together with 
the chiral extrapolations, are shown in Figure 1 
for the gauge-invariant smearing procedure. This 
procedure, with a correlated chiral extrapolation, 
provides our best value for Bg at : 
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where the systematic errors encompass the vari- 
ations due to different smearing functions, differ- 
ent fitting techniques and correlated or uncorre- 
lated chiral extrapolations. However, it must be 
emphasised that these values increase by 20 % 
if the alternative matching prescription discussed 
in Section 2.1 is adopted, and this must be re- 
garded as a further systematic uncertainty in our 
calculation. 

3. HEAVY BARYON PHYSICS 

For the study of baryons containing a single 
heavy quark, we employ propagating quarks com- 
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puted using the SW action, at four values of 
the heavy-quark mass around that of the charm 
quark, corresponding to Kh = 0.121, 0.125, 0.129 
and 0.133. Heavy quark propagators are com- 
puted from Jacobi-smeared sources to both lo- 
cal and smeared sinks; in contrast to the preced- 
ing study, we also employ light quark propagators 
from smeared sources, but only at our two heav- 
iest light-quark masses. 

3.1. Lattice operators 

There are eight baryons containing one heavy 
and two light quarks, and these are listed in Ta- 
ble 1, together with their quantum numbers and 
quark content. At lowest order in HQET it is pos- 
sible to distinguish the quantum numbers of the 
light degrees of freedom from those of the heavy 
quark, and these also are shown in the table. 

The spectrum can be computed using the fol- 
lowing operators: 

5 = e abc (l a C l5 l' b )h c (19) 
O, = Cahc(l a C^l h )h e (20) 
O'^ = e abc (l a C 7 ,l' b )h c , (21) 

where C is the charge conjugation matrix, /, /' are 
light-quark fields, and h is the heavy-quark field. 

The operator O5 corresponds to sj' = + spin- 
parity for the light degrees of freedom and a total 
spin-parity for the baryon J p = i + . We extract 
the mass from a fit to the two-point correlator 

G 5 (p,t) = Xy* 2 (0 5 (z,*)(5 5 (O,())} (22) 

X 

*>o z^_ e _ Et (23) 

The operator creates both spin ^ and spin ^ 
particles. However, projectors _P 3 / 2 and P 1 ! 2 can 
be defined [14] to separate the spin | and spin -| 
contributions to the correlator of O^, 

G,Ap,t) = ^e-^(O,(£,t)a(0,0)) (24) 

X 

2 

^L e - E i\jt + Mi)Pj v . (25) 



3.2. Masses and mass splittings 

The masses are obtained from single-state fits 
to both the LS and SS correlators of Equa- 
tions (22) and (24). In general, we find consis- 
tency between the fits, but those to the SS corre- 
lators are generally more stable under variations 
in the fitting range, and have lower x 2 /dof. We 
therefore use these for our best values. 

The heavy-baryon masses at fixed Kh are ex- 
trapolated linearly in the sum of the light-quark 
masses to the chiral limit; because we have cor- 
relators for both degenerate and non-degenerate 
light-quark masses, we can include three points 
in the extrapolation. 

Following HQET, the chirally-extrapolated 
masses are extrapolated linearly in the heavy- 
light pseudoscalar mass to obtain the various 
baryon masses at the physical D and B masses; 
the extrapolation proved insensitive to different 
modelling functions. The masses obtained from 
this procedure, together with the experimental 
values where known, are shown in Table 1. 

To focus on the dynamics of the light quarks 
and their interaction with the heavy quark, it 
is useful to study the mass splittings. Both the 
A — PS and the £ — A mass splittings are sen- 
sitive to the dynamics of the light degrees of 
freedom. Furthermore, in the latter case we ex- 
pect no dependence on the heavy-quark mass at 
0(l/mft). Below we give the dimensionless ratios 
for these splittings at both the charm and the 
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We have little confidence in our measurements 
of the spin splittings; they have large relative sta- 
tistical errors, and it is difficult to perform the 
extrapolation in the heavy-quark mass. Never- 
theless, our measurements of the spin splittings 
are generally inconsistent with experiment, and 
in most cases consistent with zero, mirroring the 
well-known feature in the meson system. 
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Table 1 



Quantum numbers of the eight baryons containing a single heavy quark h. I, sj' are the isospin and the 
spin-parity of the light degrees of freedom and S, J p are the total strangeness and the spin-parity. The 
lattice detemination is given in GeV, together with the experimental values, where known. 



3.3. Semi-leptonic form factor 

The matrix elements for the decay Aj — ► K c lvi 
are described by six form factors: 

{AUp}\J h , h '(y)\K'(f)) = 

h^\p)T h /(p',p)u { ^(f) (26) 
where 

,J h /(x) = h(x) 1 ,(l- l5 )ti(x) (27) 
and 

^• h '(p',p) = (Fii^h, + + f3("K) 

- (Gi(w) 7 ^ + G 2 (wK + G 3 (wK) 75- (28) 

Here ) are helicity indices, i/ ) are the four- 
velocities, i/'s are the heavy-baryon spinors and 
lo = v.v' . Fi(u>) and are form factors of 

the vector and axial-vector currents respectively. 
HQET requires that the form factors be related 
to a universal RGI function £(w), which we call 
the baryon Isgur-Wise function: 

F(lo) = (a,- + + 7i (a;)) (29) 

Gi(w) = (af+^H + TfHKH (30) 

a\ = a\ = 1; 0:2,3 = a| 3 = 0. (31) 

/3j 5 \ui) parametrise the radiative corrections, 
whilst jj 5 \ui) contain corrections proportional to 
powers of the inverse heavy-quark mass. 



The form factor G\{lo) has no 1/m^ corrections 
at lo = 1. Furthermore, this form factor dom- 
inates the decay rate near zero recoil, affording 
the possibility of an accurate experimental deter- 
mination of V c b- In order to compute G\{lo), we 
evaluate the three-point correlator 

(j>x i^y ) — 

^-^e-^y(Ol(x)J^ h ' (y)O h 5 ' (0)). (32) 

By taking the ratio of this to an appropri- 
ate product of two-point functions we obtain 
F^' h ' (p' , p) utt , the lattice analogue of F^' h '(p',p) 
with p 1 = p + q, and hence extract G 1 1 att (w). 

In the case of degenerate transitions, in which 
the heavy quarks are of equal mass, we can re- 
duce discretisation errors in G\{lo), as well as 
remove the dependence on the matching coeffi- 
cient Za, by dividing G 1 1 att (w) by the measured 
G\ att (l); the latter is obt ained from the transition 
(0, 0, 0) —> (0, 0, 0). In Figure 2, we show prelimi- 
nary results for G\{lo), normalised to = 1, 
at our four value of Kft , with both light quarks 
corresponding to k; = 0.14144. The error bars 
at zero recoil arise from the measurement of the 
transition (tt/12 a" 1 , 0, 0) (tt/12 a" 1 , 0, 0). 

Radiative corrections are not included in this 
preliminary analysis, so some caution has to be 
exercised before calling this an "Isgur-Wise" func- 
tion. Furthermore the power correction ji are 
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Figure 2. The form factor G\ att (u) / G\ att (1) for 
degenerate transitions, at k; = 0.14144. The dif- 
ferent plotting symbols correspond to different Kh 



not yet determined in a model-independent way. 
Nevertheless, the points appear to lie on a uni- 
versal mass-independent curve, suggesting that 
mass-dependent corrections are small; this seems 
reasonable, since G\ is free of 0(l/rrih) correc- 
tions at zero recoil. 

4. CONCLUSIONS 

The study of heavy-quark systems on the lat- 
tice is vital both as a phenomenological tool to 
enable the extraction of CKM matrix elements, 
and as a laboratory in which to explore the heavy- 
quark symmetry. In this talk, two different stud- 
ies have been presented. The measurements of 
fg at and Bg at are remarkably consistent between 
different interpolating operators and between dif- 
ferent fitting techniques. The principle require- 
ment to improve our determination still further 
is a better understanding of the matching of op- 
erators between the continuum and lattice. 

The exploratory study of the spectrum of heavy 
baryons has produced notable agreement with ex- 
periment, both in the charm and beauty sectors. 
Furthermore, we have demonstrated the feasibil- 



ity of extracting heavy-baryon semi-leptonic form 
factors, opening the possibility of further inde- 
pendent measurements of V c b- 
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